Nucleation of superconductivity in a mesoscopic loop of varying width 
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We study the evolution of the superconducting state in a perforated disk by varying the size of 
the hole. The superconducting properties are investigated by means of transport measurements 
around the superconducting/normal phase boundary Tc{H). A transition from a one-dimensional 
to a two-dimensional regime is seen when increasing the magnetic field for disks with small holes. 
A good agreement is found between the measured Tc{H) line and the calculations performed in the 
framework of the linearized Ginzburg-Landau theory. The effect of breaking the axial symmetry of 
the structure by moving the hole away from the center of the disk is also studied. An enhanced 
critical field is found for the asymmetric structures, possibly due to the recovery of the singly 
connected state. 

PACS numbers: 74.25. Dw, 74.78. Na, 74.25. Op 



I. INTRODUCTION 

The fluxoid quantization constraint in superconduct- 
ing loops structures gives rise to the well-known peri- 
odic Little-Parks^ oscillations in the Tc{H) phase line. 
The confinement of the superconducting condensate in 
mesoscopic structures leads to interesting new phenom- 
ena that are strongly dependent on the geometry and the 
topology of the structure^. 

While the limiting cases of a superconducting 
diskM^s^ ^ ^^-^ loopi^2.9.iQaia2^i3a4 ^^^^ b^^^ 

broadly studied experimentally and theoretically, the in- 
termediate case is not. The latter is also closely related 
to the problem of a thin film exposed to a parallel mag- 
netic fieldi^Ai^. In this situation, a quasi one-dimensional 
(ID) behavior, characterized by a parabolic dependence 
of the Tc{H) phase line, was predicted by Saint- James 
and de Gennes^^ as long as the thickness r (in units 
of C(T)) is smaller than approximately two. Above this 
value, the two surface superconducting sheaths are sepa- 
rated by a normal region, and a linear Tc{H) dependence 
is observed, typical for a two-dimensional (2D) system in 
a perpendicular field. At this dimensional crossover point 
(r=1.84), vortices start to nucleate in the sampleiSii^. 

In their pioneering experiment. Little and Parks^ stud- 
ied the field-temperature H — T phase diagram of a thin- 
wire loop in an axial magnetic field. A periodic compo- 
nent in the experimental Tc{H) line was found. In the 
framework of the nonlinear Ginzburg-Landau (GL) the- 
ory, Berger and Rubinstein^^'^^ studied mesoscopic su- 
perconducting loops. They have predicted that, if the 
thickness of the loop is not exactly uniform, then there 
exist situations for which superconductivity is suppressed 
at a certain location, so that the superconducting loop 
effectively becomes singly connected and no supercur- 
rent flows. When this happens, the sample is in the so 
called 'singly connected stated However, Vodolazov et al 
showed that the state where the order parameter vanishes 
at one point is still doubly connected since the phase of 



the order parameter is not independent on both sides of 
the place where the order parameter is zero. They sug- 
gested to call this state a one-dimensional vortex state. 

The intermediate case of finite width loops was studied 
within the London theory by Bardeen^^. He calculated 
that the flux is quantized in units of u^o (with u < 1) in 
cylinders of very small diameters and with a wall thick- 
ness of the order of the penetration depth. It was later 
calculated that the flux through an area S = 7rr^ is quan- 
tized in units of ^o, with r^n = (r^ +ro)/2 the arithmetic 
mean of the inner (r^) and outer (ro) radius^^. Arutunyan 
and Zharko\«Si determined in the London limit an effec- 
tive radius of Veff = ^JViVo such that inside this ring, 
the flux was exactly quantized. These two different val- 
ues and re// are nearly identical for the narrow ring. 
Baelus et alSSk found that the value of < rg// < Tq was 
dependent on the vorticity L and in fact rg// turns out 
to be an oscillating function of the magnetic field. 

A self-consistent treatment of the full nonlinear GL 
equations for a square loop has been carried out by 
Fomin et alM^. The order parameter |^| distribution 
was found to be strongly inhomogeneous due to the pres- 
ence of sharp corners. The precise shape of the Tc{H) 
curve crucially depends on the area fraction for which 
Tc{H) + 0. 

Bruyndoncx et al^^ investigated infinitely thin loops 
of finite width using the linearized GL equation. The in- 
duced magnetic field can be neglected in this case. They 
limited their investigations to circular symmetric solu- 
tions so that only the giant vortex state was studied. 
Berger and Rubinstein^^ also considered infinitely thin 
loops of finite width using the nonlinear GL theory, ne- 
glecting the induced fields. 

Baelus et al^^ analyzed circular fiat disks of nonzero 
thickness with a circular hole in it. They also investi- 
gated the case where the hole is shifted off the center of 
the disk. The superconducting properties were studied 
also deep in the superconducting state. For small super- 
conducting disks with a hole in the center, they found 



only the giant vortex state. The influence of the radius 
of the hole on the superconducting state was considered. 
For larger holes in perforated superconducting disks, a re- 
entrant behavior was seen, where a transition from the 
giant vortex state to a state with separated vortices and 
back to the giant vortex state was found. 

Recently, Pedersen et al?^ investigated experimentally 
the magnetization of a mesoscopic loop. The periodic- 
ity observed in the magnetization measurements reveals 
a sub flux quantum shift. This flne structure was inter- 
preted as a consequence of a giant vortex state nucleating 
towards either the inner or the outer side of the loop. 

This experiment has lead to a recent growing interest 
in the mechanism of flux transition in superconducting 
lQQpgi3j28^ Multiple flux jumps and irreversible behav- 
ior of the magnetization were observed in thin mesoscopic 
rings by Vodolazov et al?^ . At low magnetic fleld and 
for rings with sufficiently large radii, they showed experi- 
mentally and theoretically, using the time-dependent GL 
theory, that the vorticity may change by values larger 
than one. 

The existence of a zero-current line in mesoscopic su- 
perconducting rings has been found both theoreticall}ii& 
and through experimental observation of self-generated 
weak links^^. It was suggested in Ref. 32 that a system 
of asymmetric loops can be used as noise detector or as 
source of high frequency radiation. 

Beyond the vortex and the giant vortex conflg- 
urations, the ring-like vortex solutions of the GL 
equations in superconducting mesoscopic devices were 
investigated^SiSiS. Those solutions possess a unique 
winding number in the whole ring, but the order pa- 
rameter vanishes on one or more cylindric surface. For 
a nanosized Pb bridge, it has been reported that the 
vorticity varies along the axis of the bridge^^. Solving 
self-consistently the nonlinear GL equations for a meso- 
scopic superconducting ring, Zhao et al.^^ obtained solu- 
tions with different vorticity inside and outside the zero- 
current line at a certain radius. They, however did not 
consider the phase coupling of the order parameter be- 
tween the two superconducting parts of the ring sepa- 
rated by the zero-current line. 

The paramagnetic response for a stable configuration 
of a mesoscopic ring has been studied in Ref. They 
found an oscillation of the order parameter density profile 
when changing the vorticity. 

Using ultrasensitive susceptibility techniques and scan- 
ning Hall probe microscopy Davidovic et a/.^^'^^ have 
studied arrays of electrically isolated superconducting 
mesoscopic rings. When these rings are biased in an ex- 
ternal fiux of $0/2, they can be in either of two energeti- 
cally degenerated fiuxoid states. The magnetic moments 
produced by the supercurrents in these rings are analo- 
gous to Ising spins, and neighboring rings interact anti- 
ferromagnetically via their dipolar magnetic fields. The 
ring dynamics is dominated by an energy barrier between 
the two states which rises rapidly as the temperature is 
lowered below Tc- 



2 



(a) (b) (c) 




(d) (e) 




FIG. 1: SEM micrograph of (a) an Al disk with outer radius 
ro=l iJim and of a loop with outer radius ro=l /xm and inner 
radius (b) ri=0.1 /xm, (c) ri=0.3 /xm, (d) ri=0.5 /xm and (e) 
Ti—^.l iim. 



In this paper we shall study the systematic variation 
of the superconducting phase boundary, Tc{H)^ in per- 
forated disks with different r^/ro ratios, which realize a 
cross-over from the singly-connected disk to the limit of 
the thin ring. The rest of the paper is arranged as fol- 
lows: in Section ^ we will study the evolution of the 
superconducting state for the transition from a disk ge- 
ometry to a thin ring. The superconducting properties 
of the disks with a hole in the center will be analyzed 
by transport measurements carried out around the su- 
perconducting/normal transition line. In Section ITTll the 
effect of breaking the axial symmetry of the structure by 
shifting the hole off the center of the disk will be dis- 
cussed. The onset of dissipation below Tc{H) will be 
studied in Section HVl 



II. SUPERCONDUCTING RINGS 

A. Sample properties 

A SEM micrograph of the different studied samples 
prepared with e-beam lithography is given in Fig. ^ 
All the structures consist of disks with external radii of 
To = 1 jam. The radii of the holes (Fig. Q determined 
from SEM micrograph, were = jam (a), =0.1 /am 
(b), ri — 0.3 /im (c), = 0.5 fim (d) and = 0.7 /xm 
(e). All the samples were evaporated in the same run, 
except for the thinnest loop. A different evaporation 
will only slightly alter the superconducting properties 
like the coherence length and the critical temperature. 
Wedge shaped contacts with opening angle F = 15° are 
used in order to minimize the influence of the contacts 
on the superconducting properties of the structures^^'^^. 
The coherence length determined from a macroscopic co- 
evaporated sample was found to be ^(0) = 156 nm for 
the disk and the three loops with a small opening. The 
thickness was r = 39 nm. For the sample presented in 
Fig. de), a coherence length of ^(0) = 120 nm was de- 
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ter mined in the same way as for the other structures. A 
thickness of r = 54 nm was found from low angle X-ray 
diffraction on a co-evaporated film and from AFM for 
the loop with = 0.7 /im. The H — T phase diagram is 
constructed by four-point resistance measurements using 
an ac transport current of 0.1 /iA. The phase boundary is 
determined from a set of magnetoresistance curves mea- 
sured at various temperatures using a certain resistance 
criterion Rc. 



B. Resistance transitions 

The superconducting/normal resistance transitions for 
the disk and the rings with a different inner to outer 
radius ratio x = r^/ro are shown in Fig. [21 The five 
different samples have a very similar temperature depen- 
dence of the resistance at different magnetic field as the 
samples with wedge shaped contacts with opening angle 
r = 15° presented in Ref. ^3- They are characterized 
by a slowly decreasing resistance at high temperatures 
arising from the nucleation of superconductivity in the 
wedge contacts, followed by a sharp drop of the resis- 
tance once superconductivity nucleates in the ring. Only 
small differences are seen in the amplitude of the resis- 
tance overshoot observed at certain magnetic fields. It 
is probably due to small differences in the shape of the 
contacts that are responsible for the appearance or not 
of the resistance anomaly created by a charge imbalance 
around superconducting/normal interfaces. The samples 
with x=0.3, x=0.5 and x=0.7 show a different behavior 
at low magnetic fields. There, the situation is reversed. 
A sharp transition is firstly observed, followed by a broad 
transition at low resistance. We will show below that the 
broad transition also corresponds to the nucleation in the 
wedges. This effect is observed in a broader magnetic 
field range when the ratio x increases. 

A sharper transition at high magnetic fields is seen in 
the resistance transition if the thinnest ring. It is difficult 
to extract from the measurements if this is caused by 
the smaller coherence length or by the geometry of the 
sample. 

C. Tc{H) phase boundaries 

The experimental phase boundary of the disk is pre- 
sented in Fig. I3fa). The results are compared with 
theoretical calculations of the nucleation field H^^{T) 
(full line in Fig. by Bruyndoncx et al?^ . A very 
good agreement between the calculated and the mea- 
sured curve is seen. Only a slightly lower coherence 
length (^(0) = 130 nm) than the one found for a ref- 
erence macroscopic sample (^(0) = 156 nm) had to be 
used for the experimental data. 

The data for the ring with x = 0.1 are shown in 
Fig. Ofb). The flux ^ on the field axis denotes the 
flux $ = fioHirrl through the ring and the hole. The 



H — T diagram of the ring with the smallest hole re- 
sembles strongly the Tc{H) line of the disk displayed in 
Fig. Eta). The phase boundary has a linear background 
superimposed with oscillations. A very good agreement 
between the measured and the calculated curves is found. 

Fig. Etc) shows the H — T diagram of the ring with x = 
0.3. Here, the linear dependence is only seen for vorticity 
L > 4. At lower magnetic field, a parabolic background 
suppression of Tc is observed. The crossover from the 
linear to the parabolic regime occurs at 7rr^/^^(T) ^ 20. 
This corresponds to a value ro — Vi ^ 1.8 ^(T), which is 
in a good agreement with the thickness r = 1.84 ^(T) for 
a crossover from a ID to a 2D regime for a thin film in a 
parallel magnetic fieldiii^. 

A good agreement with the position of the cusps in the 
theoretical curve has been found. The amplitude of the 
oscillations in the experimental curve deviates slightly 
from the calculated one. At L=l, between the first and 
the second Tc{H) cusps, the experimental oscillation is 
less pronounced. For higher vorticity, the opposite sit- 
uation is seen where the amplitude of the experimental 
oscillations is larger than in the theoretical curve. 

The penetration of the first vortex in the ring occurs 
at a lower magnetic field value than for the ring with 
the smallest hole [see Fig. Ol^b)], while the transitions 
L = 1^2toL = 5^6 occur at a higher magnetic 
field. That the transitions take place at lower magnetic 
field value for a ring with thinner lines is expected since 
the transition between L and L + 1 occurs at ^/^o = 
L + 1/2 for an infinitely thin loop or cylinder. At higher 
magnetic fields, a giant vortex state is formed^ and the 
disk with a small hole in the center behaves like the disk 
without hole. This, however, cannot fully explain why 
the change in vorticity is delayed at high magnetic fields 
by introducing a small hole in a disk. 

The measured Tc{H) phase boundary of the ring with 
ratio X = 0.5 is shown in Fig. E^d). In the tempera- 
ture range accessible with our experimental setup, only a 
parabolic background dependence of the critical temper- 
ature on the magnetic field has been measured. By com- 
paring the experimental results with the calculations, a 
similar behavior as for the ring with x = 0.3 is seen. The 
position of the cusps in the experimental curve matches 
with the calculated transitions. However, no good agree- 
ment is found for the amplitude of the oscillations. For 
the vorticities L = 1 and 2, the amplitude is lower in 
the experimental curve, while for L > 3, the amplitude 
is larger. At low L, the transition between states with 
different vorticities occurs at a lower magnetic field than 
for the disk, while the transitions L = 3^4, L = 4^5 
and L = 5 ^ 6 take place at a higher magnetic field, 
similar to what was observed for the ring with x = 0.3. 

The H — T diagram of the ring with the thinnest line 
{x = 0.7) is shown in Fig.EJe). Two experimental curves 
are presented, one for Rc = 0.5 Rn (open squares) and 
the second for Rc = 0.8 Rn (open circles). It can be seen 
that at a higher resistance criterion the parabolic depen- 
dence switches to a linear regime at high magnetic field. 
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FIG. 2: Resistive transitions R{T) for (a) a disk and for a loop with inner to outer radius ratio (b) x = 0.1, (c) x = 0.3, (d) 
X = 0.5 and (e) x = 0.7 in different magnetic fields. The dashed and dotted lines show the resistance criteria used to determine 
the Tc{H) phase boundary. 



For the curve calculated with the low resistance criterion, 
a quasi-parabolic background suppression of Tc{H) is ob- 
served over the whole measured range. The amplitude of 
the Tc{H) oscillations is larger than in the samples with 
smaller x and the transition between states with different 
vorticities is almost periodic in field. A good agreement 
between the theoretical curve and the experimental curve 
with Rc = 0.5 Rn is seen at high magnetic fields. At lower 
magnetic fields, a good agreement is found when using a 
higher resistance criterion. 

The phase boundaries of the four different loops are 
compared with the critical temperature of the disk in 
Fig. Ol^f). All the curves overlap with each other for 
L = 0. R is interesting to note that an opening in the disk 
does not affect the phase boundary as long as no vortex 
is trapped inside the superconductor. Only the magnetic 
field range over which the state with L = exists at the 
phase boundary is lowered by introducing a hole in the 
disk. The Tc{H) line of the disk with the smallest hole in 
the middle does not deviate substantially from the phase 
boundary of the disk without any opening. Only small 
changes in the position of the cusps is observed at low 
vorticity. For larger holes, the crossover from 2D to ID 
regime is clearly seen. The samples with the thinnest 



lines do not show the 2D regime in the studied tempera- 
ture interval and only the parabolic dependence is seen. 

In order to reveal the origin of this different behav- 
ior at low and high magnetic field, a contour plot of the 
resistance R{H^ T) is presented in Fig. 01 Two differ- 
ent parts are clearly distinguished. Below 10 mT, the 
low resistance is linear, while the high resistance exhibits 
a parabolic background superimposed with oscillations. 
Above 10 mT, the opposite situation occurs, where the 
low resistance has a parabolic decay with small oscilla- 
tions while the high resistance decreases monotonously. 
The parabolic part coincides with the nucleation of super- 
conductivity in the loop shown as a full line. The linear 
part arises from the nucleation in the wedge contacts. 

By fitting the theoretical critical temperature of a 
wedge with opening angle F = 15° to the linear part 
of the contour plot (dashed line), a coherence length 
^(0) = 140 nm is obtained. This differs from the coher- 
ence length ^(0) = 110 nm that was used to find a good 
agreement between the experiment and the theoretical 
curve of a loop. A possible origin of this discrepancy 
could be a width of the loop that has been evaluated to 
be smaller than the real size. An estimate of the thick- 
ness that would satisfy the coherence length used for the 



5 




30 
25 
20 
15 
10 
5 



(d) 



R =0.5 
a(Q)=139nm 
r.=1.01 |jm 




1 23456789 10 





6 8 10 12 14 16 18 



FIG. 3: Experimental Tc{H) phase boundary of (a) a disk and for a loop with inner to outer radius ratio (b) x = 0.1, (c) 
X = 0.3, (d) X = 0.5 and (e) x = 0.7. The open squares [and the open circles in (e)] represent the measured data. The 
experimental Tc{H) phase boundaries of the different structures are compared in (f). 
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FIG. 4: Contour plot of the resistance R{H,T) of a loop 
with X = 0.7. The full line represents the calculated phase 
boundary of a loop with n = 0.7 /urn and To = 1 /xm, using a 
coherence length of 110 nm. The dashed line is the theoretical 
critical temperature of a wedge with opening angle F = 15° 
with ^ 140 nm. 

calculation of the wedge contacts can be obtained from 
the analysis of the nucleation field of a thin wire of a 
film in a parallel magnetic field. From the calculation of 



the nucleation field of a thin film in a parallel field^^^, a 
value for the width of the loop of 0.38 /im is obtained 
instead of 0.3 /im found from SEM measurements. This 
difference is too large to be explained only by an error 
in the characterization of the sample. The opening angle 
of the contacts can be determined with a high accuracy 
so that a divergence arising from a wrong determina- 
tion of F could be excluded. It means that either the 
nucleation of superconductivity is delayed in the wedges 
due to the presence of the loop or that the nucleation in 
the loop is enhanced by the contacts. It is also possible 
that the coherence length in the loop is slightly different 
from that in the wedge. The sample geometry can in- 
deed affect the superconducting parameters A and ^ in 
a structure of mesoscopic size similar to the case of a 
thin film where the effective penetration depth increases 
as A' = A^/r, taking into account the demagnetization 
effects. The renormalization of A and ^ should therefore 
be calculated in a self-consistent way from the sample 
geometry. 

The shape of the resistive curves in Fig.I2fe) can be eas- 
ily understood from Fig.^ It was clearly seen that in low 
magnetic fields the nucleation first occurs in the ring and 
is then followed by the nucleation in the contacts. Due 
to the different field dependence of the Tc{H) of the ring 
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and the contacts, the opposite occurs in higher magnetic 
fields. Two different shapes are therefore distinguished 
in the resistive curves depending on the part where su- 
perconductivity starts to nucleate. The same happens in 
the rings with x = 0.3 and x = 0.5 since Tc also has a 
parabolic field dependence for low fields so that the broad 
transition at low resistance seen in Figs.EIJc) andElJd) is 
due to the nucleation of superconductivity in the wedge 
that takes place after the nucleation in the ring at low 
magnetic fields. The normal parts of the sample can how- 
ever partially become superconducting by the proximity 
effect with the neighboring superconducting part. 



D. Resistance criterion 

Since the contacts have a different field dependence 
than the studied structures, the choice of a constant resis- 
tance criterion for the determination of the Tc{H) phase 
boundary seems not to be obvious anymore. In order to 
study the best resistance criterion that should be used, 
the experimental resistance R{H^ T) was determined at 
the temperature T{H) corresponding to the theoretically 
calculated phase boundary Tc{H). The results are given 
in Fig. [51 for the loops with x = 0.3, x = 0.5 and x = 0.7. 
The curve corresponding to the resistance of the ring with 
the thinnest hole {x = 0.1) is not shown. These curves 
show the expected behavior with an approximately con- 
stant resistance equal to the resistance criterion used for 
the determination of Tc{H) presented in Fig. IHl^b). For 
other rings a similar R{H, T) dependence was obtained. 
At low magnetic field, the resistance is high, and it drops 
to a low value above a certain magnetic field. The field 
where this transition occurs increases with increasing x 
and corresponds approximately to the position where the 
phase boundary of the contacts crosses the phase bound- 
ary of the loop. For the ratio x = 0.3, a normal regime 
is found at high magnetic field with an almost constant 
resistance with a value approximately equal to the resis- 
tance criterion used for the determination of the phase 
boundary. It occurs when the linear regime is recovered. 
This is not seen in the two other samples since there, 
the linear regime is not attained. At low magnetic field, 
superconductivity nucleates first in the ring. This is the 
upper part in the resistance curves. At this point, a high 
resistance criterion should be taken. For higher magnetic 
fields, the resistance starts to drop once the contacts be- 
come superconducting. In order to determine the H — T 
diagram of the ring in that region, a low resistance tran- 
sition should be taken. This is exactly what is seen in 
Fig.O The exact shape of the curves shown in the figures 
strongly depends on the coherence length used for the 
calculation of the theoretical Tc{H) line, but the general 
behavior will not strongly change while using a slightly 
different value of ^(T). 

The origin of the discrepancy in the amplitude of os- 
cillations in the Tc{H) phase boundary is most probably 
due to the fact that a constant criterion works well to 



determine the phase boundary when the critical temper- 
ature of the contacts has a similar field dependence than 
the studied sample geometry. When this is not the case, 
the determination of the phase boundary is strongly hin- 
dered. 

It is also interesting to note that oscillations are present 
in the curves, with a maximum at the transition between 
different vorticities which becomes a minimum at high 
magnetic fields. This crossover also corresponds approx- 
imately to the field where the Tc{H) curves of the loops 
and of the contacts cross each other. A minimum in the 
resistance curves of Fig. El is observed when the differ- 
ence in critical temperature of the loop and of the wedge 
contacts is minimal and a maximum when the difference 
is maximal. This also refiects the observed difference in 
amplitude of the oscillations in the experimental and the- 
oretical phase boundaries. 

Meyers^^ calculated that the order parameter evolves 
from stronger at the inner part to stronger at the outer 
part when the vorticity is increased. This could affect our 
measurements since our applied transport current would 
fiow along the inner edge at fields slightly lower than 
the field where the transition of vorticity happens and 
then fiow along the outer edge for slightly higher fields. 
Since our contacts are situated on the outer edge, a higher 
resistance could be expected when the transport current 
is fiowing along the inner edge. We however do not see 
an increase (decrease) of the resistance just before (after) 
the cusp. 



III. NON-SYMMETRIC GEOMETRIES 

In this section the nucleation of superconductivity is 
studied for disks with a hole. The aim of this study is to 
analyze the effect of displacing the hole from the middle 
of the disk on the phase boundary. 

A SEM micrograph of the studied structures is pre- 
sented in Fig. El The three samples are the disks with 
outer radius To = 1 /im and contain a circular open- 
ing with radius = 0.3 /im. The hole is in the center 
of the disk [Fig. Eta)] or is displaced from the center 
of the disk over a distance a = 0.3 /im [Fig. Etb)] and 
a = 0.6 jiim [Fig. Etc)]. The two asymmetric samples 
were co-evaporated in the same run as the circular sym- 
metric ring with ratio x = 0.3 discussed in Section 
They have the same thickness r = 39 nm. The coher- 
ence length of ^(0) = 156 nm was determined from a 
reference macroscopic film. Wedge shaped current and 
voltage contacts with an opening angle F = 15° were 
used. 

The resistive transitions of the rings with a = 0.3 /xm 
and a = 0.6 fim are shown in Fig. [3 The transitions at 
2 mT for a = 0.3 /im and at 1 and 2 mT for a = 0.6 /im 
exhibit a sharp drop by decreasing the temperature fol- 
lowed by a slowly decaying resistance in the lowest part 
of the curves. This is similar to the curves of the sym- 
metric ring shown in Fig. I2tc), but less pronounced. It 
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FIG. 5: Resistance R{H,T) of a loop with (a) x — 0.3, (b) x = 0.5 and (c) x = 0.7 measured at a temperature T{H) = Tc{H) 
following the theoretical phase boundary Tc{H) shown in Fig. EIc), (d), and (e), respectively. The dotted lines represent the 
position of the cusps in the theoretical H — T diagram. 
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FIG. 6: SEM micrograph of an Al loop with outer radius 
ro=l /im and inner radius n^O.S /am with the hole (a) in the 
center, (b) moved over a distance a — 0.3 jim and (c) over 
a = 0.6 iim from the center. 



was seen in Section ^ that the sharp decrease of R is due 
to the nucleation of superconductivity that occurs first 
in the ring at these magnetic fields. For higher fields, 
the curves show transitions that are similar to the R{T) 
curves of the disk. The sharp part at /io^ ^ 3 mT also 
corresponds to the nucleation of superconductivity in the 
ring. At higher magnetic fields, the nucleation starts first 
in the wedges and is then followed by the ring. 

The H — T diagram of the disk with an off-centered 
hole displaced by a = 0.3 iim off the center is shown 
in Fig. [Sfa). Small oscillations are seen in the phase 
line. A behavior in between the parabolic and the linear 
field dependence is observed. At higher magnetic field, 
the oscillations are almost not distinguishable anymore. 
The phase boundary of the disk with the hole displaced 
by a = 0.6 jim from the center is given in Fig. IHI^b). 
In the temperature range accessible in our experimental 
setup, only a linear regime was observed. Also very weak 
oscillations were distinguished. 

The two Tc{H) lines are compared with the phase 
boundary of the circular symmetric ring, with the same 
inner and outer radii, shown in Fig. IHIJc). The three 
curves have approximately the same behavior at low 
magnetic field, except that the oscillations are less pro- 



nounced for the non- symmetric rings. It was seen in Sec- 
tion^ that an opening does not affect the phase bound- 
ary for L = 0. There, the symmetry was kept. When 
the circular symmetry is broken, the phase boundary is 
strongly affected as can be seen in the inset of Fig. [SJc) 
even when no vortices are trapped in the sample. At 
high magnetic fields, the curves separate. Increasing the 
asymmetry enhances superconductivity. Baelus et al?'^^ 
calculated the free energy, the magnetization and the 
Cooper-pair density of non-symmetric rings with finite 
width. They found that the density of the superconduct- 
ing condensate was the highest in the narrowest region of 
the superconductor when L 7^ 0. They argued that the 
trapped flux tries to restore the broken symmetry. That 
superconductivity is stronger in the smallest part of the 
sample is probably due to the fact that the critical field 
is enhanced in thin lines. With the configuration of the 
contacts that was used (see Fig. a superconducting 
'bridge' can be formed across which the external current 
applied for transport measurements can pass. The criti- 
cal field of this area will probably be higher than in the 
one with the largest area of superconducting material. 
The measured phase boundary is therefore most likely 
only the phase boundary of the bridge and not of the 
full sample. Given that no supercurrent can circulate 
around the opening, a singly connected state is then re- 
covered^^ . Since a superconducting path will always 
be found across the bridge, a lower resistance criterion 
for Tc{H) will not determine the nucleation in the whole 
sample. The phase boundary of the complete structure 
could only be probed with contacts turned by 90°. 

It is worth emphasizing that the phase boundary of the 
structure with a = 0.6 /im exhibits a linear field depen- 
dence, typical for a 2D behavior. The thinnest part of the 
sample where the nucleation firstly occurs can be seen as 
a curved line of varying width. When the asymmetry is 
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FIG. 7: Resistive transitions R{T) in different magnetic fields for a ring with outer radius ro=l /xm and with a hole radius 
ri=0.3 /jm. The hole is moved (a) by 0.3 /jm and (b) by 0.6 /jm from the center. The dashed line shows the resistance criterion 
used to determine the Tc{H) phase boundary. 







FIG. 8: Experimental Tc{H) phase boundary of a loop with outer radius ro=l /xm and with a hole radius r^^O.S /am determined 
for a resistance criterion of 1/2 Rn- The hole is moved (a) by 0.3 /am and (b) by 0.6 /xm from the center. The open squares 
represent the data normalized by the coherence length ^(0) = 135 nm. (c) Comparison between the different phase boundaries. 
The hole is in the center (open squares) or moved by 0.3 /xm (open circles) and 0.6 /am (open triangles) from the center. The 
inset shows a magnification of the low magnetic field region. 



less pronounced (a = 0.3 /im), a curve in between the 
parabolic and the linear regime is observed. This indi- 
cates that the superconducting path resembles a thin line 
but can not be fully considered as such. 

The transition from L = ^ 1 is strongly delayed 
when increasing the asymmetry. The first vortex enters 
the sample at ^ = 1.3 $o for the symmetric sample and 
at <l> = 1.4 ^0 and ^ = 2.9 for a = 0.3 and 0.6 /im 
respectively. Baelus et a/.^^ indeed found a delay of the 
penetration of the first vortex by increasing the displace- 
ment of the hole from the center but not with a factor 
larger than two as seen for the sample with a = 0.6 /im. 
The next transition (L = 1 ^ 2) is not clearly resolved. 
We could however see that the transition occurs first for 
the sample with a = 0.3 /im at $ 3.4 <l>o, followed 
by the symmetric ring at ^ = 3.7 <l>o, and then only 
at <l> = 4.7 <l>o for the strongly asymmetric sample, i.e. 
the sample with the largest off the center shift of the 



hole. The theoretical investigations, however, showed a 
decreasing field for the entry of the second vortex by in- 
creasing the asymmetry. The calculations were not per- 
formed at the phase boundary but deep in the super- 
conducting state. Also a slightly smaller hole than in our 
experiments and different material parameters were used. 
This could be a reason for the discrepancy for the sample 
with a = 0.6 /im. Deep in the superconducting state, the 
complete structure will be in the superconducting state, 
while at the measured phase boundary, only the region 
around the hole is superconducting so that a smaller ef- 
fective area should be taken, to explain the strong in- 
crease of the magnetic field value where the transitions 
between L ^ L -\- 1 take place in the sample with the 
hole displaced over the largest distance from the center. 

Baelus et a/.^^ found that at L = 1 the vortex was 
trapped in the hole and that at L = 2 and 3, the hole 
captures one <l>o-vortex. A second vortex (2$o-vortex for 



9 



1/ = 3) is placed across the axis of displacement of the 
hole but at the opposite position. In our experiment, we 
believe that the sample area where the second vortex was 
found in the calculation is not yet in the superconducting 
state. In this case, no vortex can sit there and the hole 
can also not trap any vortex since no superconducting 
path around the hole exists. The vortices must then be 
placed at a position between the contacts and the hole. 
This configuration is similar to the case of a infinite wedge 
where confined circulating supercurrents were predicted 
in the vicinity of the corner^^. 



IV. DISSIPATION BELOW Tc{H) 

In order to analyze the onset of dissipation below 
Tc{H), the resistance has been measured at a certain 
fixed temperature below the phase boundary. Assuming 
that the resistance criterion {Rc = 2/3Rn) used for the 
determination of the phase diagram of the disk presented 
in Fig. Ol^a) is correct, the experimental phase bound- 
ary was shifted along the temperature axis with a value 
AT and the resistance has been measured following the 
translated phase boundary, i.e. at a temperature equal to 
T{H) = Tc{H) - AT. These curves were obtained from 
the same set of R{H) curves as for the determination of 
the phase boundary. The result for the disk is shown in 
Fig. Eta). The curve with zero shift gives the resistance 
criterion of 2/3 as expected. For temperatures below 
the phase boundary, two different parts are directly dis- 
tinguished. For low fields, an almost zero resistance is 
measured, even for the smallest shift of the phase bound- 
ary. At higher fields, a resistive region is found. The re- 
gion where no resistance is observed corresponds to the 
magnetic field value where the vorticity is zero in the 
disk {fioH < 1.3 mT). Once that the first vortex enters 
the sample a resistive behavior is noticed. The reason 
for the observation of a finite resistance can be found in 
the presence of the external current used to measure the 
phase boundary. Just below the critical temperature of 
the disk, the value of the order parameter is very low so 
that a low current will easily destroy locally superconduc- 
tivity. When going deeper in the superconducting state, 
a larger area can sustain the applied current up to a cer- 
tain temperature where a superconducting path over the 
sample is found and a zero resistance is measured. This 
however can not fully explain the magnetic field depen- 
dence shown in Fig. Efa). Another possible dissipation 
mechanism, also related to the applied current, could be 
the motion of the vortices. A current / will generate a 
Lorentz force T^ ex / x <5o on the vortices. The electric 
fields generated by the vortex motion can cause dissipa- 
tion of energy that is characterized by the observation of 
a finite resistance in transport measurements. The dissi- 
pation can also be caused by the nucleation of phase slips 
centers in the sample. It was found in^^ that the change 
of vorticity in a superconducting loop transit through a 
phase slip state associated with a smaller energy barrier 



for the transition from T to T ± 1. 

Pronounced oscillations are present in the resistance 
curves of the disk. While in magnetoresistance curves 
at a constant temperature, the appearance of oscillations 
is directly related to the presence of cusps in the Tc{H) 
line. Here no oscillations should be expected since the 
resistance is measured at a fixed temperature interval 
below the phase boundary. The position of the peaks 
in Fig. Ol^a) corresponds to the magnetic fields where 
cusps are observed in the phase boundary [dotted lines in 
Figs.Ol^a) andini^a)]. The existence of the oscillations sug- 
gests that the dissipation mechanism strongly depends on 
the stability of the vortices. The large amplitude of the 
oscillation shows that it is easy to move the vortices only 
when the magnetic field is close to the value where the 
vorticity changes: L ^ L -\- 1. At these field values, a 
constant fiuctuation between the states with vorticity L 
and T + 1 will probably occur. Vortices will then enter 
from one side and will leave the sample at the other side 
in the direction imposed by the Lorentz force. At high 
vorticity a large dissipation is observed indicating that 
the motion of vortices is more pronounced when more 
vortices are present in the sample. 

The resistance of the rings measured at a tempera- 
ture T{H) = Tc{H) — AT below the experimental phase 
boundary is shown in Fig [U| for the rings with a = 0, 0.3 
and 0.6 /im, respectively. The dissipation in the circular 
symmetric ring resembles that of the disk. Peaks are seen 
at the magnetic field value where the transition between 
two states with different vorticity takes place. The resis- 
tance at a fixed temperature interval below the critical 
temperature seems to be smaller than for the disk with- 
out opening. The hole in the center of the disk will most 
probably act as an artificial pinning center, preventing 
the vortices to move. The dissipation in the asymmetric 
samples starts to grow as soon as one vortex enters the 
sample, as in the symmetric sample. However, for larger 
magnetic fields, the two asymmetric samples have a very 
different behavior. The oscillations seen in the curves of 
the symmetric sample are almost completely suppressed. 
Moreover, the resistance is not continuously growing but 
seems to saturate above a certain magnetic field. For 
AT = 20 mK, the resistance even decreases with increas- 
ing magnetic fields for the sample with a = 0.3 jam. The 
dissipation is much lower than for the symmetric sample. 
All these observations indicate the presence of a quite dif- 
ferent mechanism of the vortex motion in the two asym- 
metric samples compared to the symmetric structures. 
Unfortunately, all the measurements were performed us- 
ing an ac current. By comparing the dissipation mea- 
sured with dc current in two directions, it could be pos- 
sible to detect if a preferential trajectory for the vor- 
tex motion exists. Since the samples are not symmetric 
around the line between the two current contacts, the 
Bean-Livingston barrier should also be asymmetric. It is 
then natural to expect a vortex motion that is dependent 
on the sign of the applied current, so called 'vortex diode' 
effect^ 
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FIG. 9: Resistance of (a) the disk and of the ring with outer radius ro=l fJ^m and with a hole radius ri=0.3 fim. The hole is 
(b) in the center and moved (c) by 0.3 iim and (d) by 0.6 i^m from the center. The resistance is measured at a temperature 
T{H) = Tc{H) — AT corresponding to a shift of AT = 0, 5, 10, and 20 mK below the experimental phase boundaries. 



V. CONCLUSIONS 

We have studied the nucleation of superconductivity 
in doubly connected superconductors in the form of thin 
superconducting disks with a circular opening. The ef- 
fect of the size and of the position of the hole on the 
superconducting properties of the structures has been in- 
vestigated. A parabolic background of Tc with periodic 
oscillations is found for the thinnest loops. For disks 
with smaller holes, a transition from a ID regime to a 
2D regime is seen when increasing the magnetic field. 
For high magnetic field, the loops recover the behavior 
of the disk without opening. A giant vortex state is then 
formed and the opening in the middle of the disk does 
not play an important role anymore. 

The experimental results of the rings of different wire 
width were compared with theoretical calculations in the 
framework of the linearized GL equation. Good agree- 
ment between our experimental results and the calcula- 



tion of Tc{H) were found. Small deviations in the ampli- 
tude of the oscillations were observed. Moreover, for the 
thinnest loop that was studied, two different resistance 
criteria had to be used for low and high magnetic fields. 
These deviations were explained by the fact that the the- 
oretical linear phase boundary of the contacts is crossing 
the parabolic Tc{H) line of the thin loops. At low mag- 
netic fields, the nucleation occurs first in the loop while 
at higher magnetic fields superconductivity develops first 
in the wedge shaped current and voltage contacts. As 
a consequence, a resistance criterion, dependent on the 
magnetic field, should be used for the determination of 
Tc{H). 

Breaking the symmetry by moving the hole away from 
the center increases the critical field. The displaced 
hole forms a small region where superconductivity is en- 
hanced. A superconducting path for the applied cur- 
rent is likely to be formed before superconductivity nu- 
cleates in the whole sample. The super current cannot 
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flow around the hole so that the singly connected state 
is recovered for a loop. The dissipation mechanism due 
to vortex motion is strongly altered in this case. 

It has been observed that the phase boundary is not 
affected by the presence and by the dimensions of a hole 
as long as no vortex is trapped inside the sample. This, 
however, is only valid when the circular symmetry of the 
structure is kept. Once that the symmetry is broken by 
shifting the hole from the center, the phase boundary at 
L = deviates from the Tc{H) line of a disk. 
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